Non-Markovian band-edge effect and entanglement generation of quantum dot 
excitons coupled to nanowire surface plasmons 
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The radiative decay of quantum dot (QD) excitons into surface plasmons in a cylindrical nanowire 
is investigated theoretically. Maxwell's equations with appropriate boundary conditions are solved 
numerically to obtain the dispersion relations of surface plasmons. The radiative decay rate of QD 
excitons is found to be greatly enhanced at certain values of the exciton bandgap. Analogous to 
the decay of a two-level atom in the photonic crystal, we first point out that such an enhanced 
phenomenon allows one to examine the non-Markovian dynamics of the QD exciton. Besides, due 
to the one dimensional propagating feature of nanowire surface plasmons, remote entangled states 
can be generated via super-radiance and may be useful in future quantum information processing. 
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The collective motions of an electron gas in a metal 
or semiconductor are known as plasma oscillations. In 
the presence of surfaces, not only the bulk modes are 
modified, but also surface modes can be created 1 . When 
a light wave strikes a metal surface, a surface plasmon 
polariton — a surface electromagnetic wave that is cou- 
pled to plasma oscillations, can be excited. Investi- 
gations of the dispersion relations of surface plasmons 
for different geometries have been reported 2 since the 
1970s. Recently, great attention has been focused on 
the so-called plasmonics since surface plasmons reveal 
strong analogies to light propagation in conventional di- 
electric components 3 . For examples, it is now possible 
to confine them to subwavelength dimensions 4 and de- 
velop novel approaches for waveguiding below the diffrac- 
tion limit 5 . The useful subwavelength confinement, sin- 
gle mode operation 6 , and relatively low power propaga- 
tion loss 7 of surface plasmon polaritons could be applied 
to miniaturize photonic circuits 8 . High surface plasmon 
field confinement was also used to demonstrate an all- 
optical modulator 9 . 

Plasmon induced modification of the spontaneous 
emission (SE) rate is naturally an extended issue. 
Arnoldus et al. theoretically investigated the atomic flu- 
orescence and spontaneous decay near a metal surface 10 . 
Strong enhancement of fluorescence due to surface plas- 
mons was also observed 11 . This enhanced fluorescence 
was considered as a possible method to improve the quan- 
tum efficiency of light-emitting devices. R. Paiella re- 
cently proposed tunable surface plasmons in silver-GaN 
multiple layers to increase the radiative recombination 
rate and equivalently the photoluminescence efficiency 12 . 
Strong and coherent coupling between individual optical 
emitters and guided plasmon excitations in conducting 
nanowires at optical frequencies was also pointed out 13 
and may be used as a novel single-photon transistor 14 . 



In this work, we investigate the SE rate of a II- VI col- 
loidal QD (nanocrystals) exciton coupled to surface plas- 
mons in a silver nanowire. Radiative decay of a QD exci- 
ton into different modes of surface plasmons is considered 
separately. The emission rate is found to be greatly en- 
hanced at certain values of QD exciton bandgap, which 
is similar to the band-edge effect in photonic crystals. In 
addition, application of such a system in generating re- 
mote entangled states is also pointed and may be useful 
in future quantum information processing. 



Dispersion relations of surface plasmons. -Consider 
now a colloidal CdSe/ZnS QD near a cylindrical silver 
nanowire with radius a. The QD and nanowire are as- 
sumed to be separated by a GaN layer 15 as shown in Fig. 
1 (a). One of the main reasons to choose a CdSe/ZnS QD 
exciton as the two-level emitter is that it is now possi- 
ble to isolate single colloidal QD and measure its exciton 
lifetime 16 . The other reason is that its exciton bandgap 
is around 2eV to 2.5eV, depending on the size and envi- 
ronment of the dot 17 . The plasmon energy hu> p of bulk 
silver is 3.76 eV with the corresponding saturation en- 
ergy hujp/^/2 s» 2.66eV in the dispersion relation 18 . As 
we shall see below, variations of the dispersion relations 
in energy just match the exciton bandgap of colloidal 
CdSe/ZnS QDs. 



Surface plasmon modes arc created due to the nonzero 
local charge density on the surface of a nanowire. The 
n-th surface plasmon mode's components of the electro- 
magnetic field at the surface can be obtained by solv- 
ing Maxwell's equations in a cylindrical geometry (p and 
tp denote the radial and azimuthal coordinates, respec- 
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FIG. 1: (a) Schematic view of the model: a silver nano-wire 
is embedded inside GaN material and a QD (quantum dot) is 
put on top of it. (b) Dispersion relations of surface plas- 
mons for the first few modes. The units for vertical and 
horizontal lines are Q = uj/u) v and K = k z c/u p . (c) Cor- 
responding (Markovian) SE rates into surface plasmons. (d) 
Non-Markovian decay dynamics of QD excitons for 5 — 0.2/3 
(black line), 0.4/3 (dotted line), and 0.8/3 (dashed line). As 
5 — 0, the red line represents the result for the contribution 
from n — 1 mode. 



tively) with the appropriate boundary conditions 2 : 
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H z = [^i(K iP )Bi}cb n , 
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where J n {Kj P ) are H„ (Kop) are Bessel and Hankel 
functions, respectively. I (O) stands for the component 



inside (outside) the wire. The dielectric function is as- 

sumed as e(u>) = Eqo[1 — u(u+iM h wnere £ oo = 9.6 (for 
Ag), Coo = 5.3 (for GaN), and r is the relaxation time due 
to ohmic metal loss 12 . The magnetic permeabilities Pl 
are unity everywhere since we consider nonmagnetic ma- 
terials here. A£ and B^ are constants to be determined 
by normalizing the electromagnetic field to the vacuum 
fluctuation energy, J e(\E p \ 2 + \E^\ 2 + \E z \ 2 )dr = huj(k), 
and matching the boundary conditions. The dispersion 
relations of the surface plasmons are thus obtained by 
solving the following transcendental equation numeri- 
cally: 
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The dispersion relations for various modes n are shown 
in Fig. 1 (b) with effective radii R = 0.1. One unit of the 
effective radii R (= uj p a/c) is roughly equal to 53.8 nm. 
The behavior for the n — mode is very similar to the 
two-dimensional case 18 , i.e. Q, gradually saturates with 
increasing wave vector k z . This is because the fields for 
the n — mode are independent of the azimuthal angle 
tp. However, the behavior for the n ^ modes are quite 
different. The first interesting point are the discontinu- 
ities around lu/c w k z . Further analysis shows that the 
solutions of oj are "almost real" 19 as k z > Re[u>]/c. Thus, 
the first Hankel function of order n, Hn\K^p), decays 
exponentially. This means the surface plasmons in this 
regime are confined on the surface (bound modes). For 
k z < Re[uj]/c, however, the solutions of uj are complex 
as shown by the dashed lines in the inset of Fig. 1 (b). 
H]i (K^p) in this case is like a traveling wave with finite 
lifetime (non-bound modes). 

Once the electromagnetic fields are determined, the 
spontaneous emission (SE) rate, r sp , of the QD exci- 
tons into bound surface plasmons can be obtained via 
Fermi's golden rule. The SE rates of the first few modes 
(n = 0, 1,2, 3) are shown in Fig. 1 (c) with effective radii 
R = 0.1. In plotting the figures, the distance between 
the dot and the wire surface is fixed as d = 10.76 nm. 
The novel feature is that the SE rate approaches infin- 
ity at certain values of the exciton bandgap uiq. Mathe- 
matically, one might think that at these values the corre- 
sponding slopes of the dispersion relation are zero. Phys- 
ically, however, this infinite rate is not reasonable since 
it's based on perturbation theory. Therefore, one has 
to treat the dynamics of the exciton around these val- 
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ues more carefully, i.e. the Markovian SE rate is not 
enough. One has to consider the non- Markovian be- 
havior around the band-edge, which means the band 
abruptly appears/disappears across certain values of u>. 

Non- Markovian dynamics of QD excitons. -To obtain 
the non-Markovian dynamics of the exciton, we first write 
down the Hamiltonian of the system in the interaction 
picture (with the rotating wave approximation), 
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= \i) =1,1) are the atomic operators; 

a n ,kz and a n k are the radiation field (surface plasmon) 
annihilation and creation operators; An,k, = ^n.k z ~ w o 
is the detuning of the radiation mode frequency 0J n ,k z 
from the excitonic resonant frequency loq, and g n ,k z = 

A 1 ■ £n,t, is the atomic field coupling. Here, ~jl and E n ^ z 
denote the transition dipole moment of the exciton and 
the electric field, respectively. 

Assuming there is an exciton in the dot with no plas- 
mon excitation in the wire initially, the wavefunction of 
the system then has the form 



iv(t)> = &e(*)ir,o>+x; w*)u>i* 
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The state vector |"f , 0) describes an exciton in the dot and 
no plasmons present, whereas ||, l n ,0 describes the ex- 
citon recombination and a surface plasmon emitted into 
mode k z . With the time-dependent Schrodinger equa- 
tion, the solution of the coefficient b e (t) in z-space is 
straightforwardly given by 
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In principle, b e (t) can be obtained by performing a nu- 
merical inverse Laplace Transformation to Eq. (5). 

To grasp the main physics and without loss of gen- 
erality, we focus on the values of u>q close to one par- 
ticular local extremum. In this case, the dispersion rela- 
tion for this particular n mode around the extremum can 
be approximated as LUk z — Lo n _ c ± A n (k z — fc„. c ) 2 ,where 
the extremum is located at (A;„ jC , (j„ i( .), The +/— sign 
represents the approximate curve for the local mini- 
mum/maximum of the dispersion relation. Once we make 
such an approximation, the radiative dynamics of the QD 
exciton is just like that of a two-level atom in a photonic 
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where S = ujq — uj ntC is the detuning and 7 is the decay 
rate contributed from other modes. 

The coefficient b e (t) can now be obtained by perform- 
ing the Laplace transformation to Eq. (5) 20 . The black, 
dotted, and dashed lines in Fig. 1 (d) represent the de- 
cay dynamics of the QD excitons for different detunings: 
S = 0.2(3, 0.4/?, 0.8/3, respectively. Here, (3 is the decay 
rate of the QD exciton in free space. In plotting the fig- 
ure, ujq is chosen to be close to the local minimum of the 
dispersion relation of the n = 1 mode. The radius of the 
wire and the wire-dot separation are identical to those 
in Fig. 1 (b). As can be seen, there exists oscillatory 
behavior in the decay profile, demonstrating that decay 
dynamics around the local extremums is non-Markovian. 
If one considers only the contribution from the n = 1 
mode and set the detuning 5 — 0, the probability ampli- 
tude would saturate to a steady limit as show by the red 
line. This quasi-dressed state is an analogous of Rabi- 
oscillation in cavity quantum electrodynamics, and also 
appears in the systems of photonic crystals 20 . 
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FIG. 2: Variations of Re[b 2 (t)] (solid lines), Im[& 2 (t)] (dashed 
lines), and |&2(i)| 2 (inset) as functions of time for ujq = 
0.602oj p (red lines) and 0.748oj p (black lines). In plotting the 
figure, the inter-dot distance zo is set equal to 0.35(w p a/c) 
with radius R = 0.1. 

Application in entanglement generation. -Let us now 
put another QD close to the wire, the interaction between 
the wire and QDs can now be written as 

H ' = £ 9n,k z (0-1+ + o- 2 +e ik * Zo )a njkz + h.c, (7) 

n : k z 

where a n ,fe s is the surface-plasmon operator, Uj + is the 
creation operator of the j-th QD. Note that in Eq. (7) 
we have assumed no detuning and the two dots have the 
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same separation from the metal wire. Since the propa- 
gating modes are along the z-direction only, the phase 
difference acquired by the second dot is ik z z , where z 
is the separation between the two dots. If one further as- 
sumes that only QD-1 is initially excited, the state vector 
of the system can be written as 
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with 61 (0) = 1 and 6 2 (0) = 6„, fea (0) = 0. Here, |||;0) 
(||t; 0)) means that QD-l(-2) is excited, while \H; l n ,fcj 
represents that both the QDs are deexcited with the pres- 
ence of single surface-plasmon. If we let the exciton band- 
gap loq far away from the band-edge, b\{t) and b%(t) 
can be obtained easily by solving the time-dependent 
Schrodinger equation. 

Fig. 2 shows the time variations of Ke[b2(t)}, lm[b2(t)], 
and |fo2(i)| for different values of the exciton bandgap. 
For wo = 0.748w p , the population of the second dot 
vanishes quickly as seen from the black lines. For 
luq = 0.602oj p , however, it approaches (quasi-)stationary 
limit 21 as shown by the red lines. This is because, for 
the later case, only n = mode contributes to the decay 
rate T sp . In this case, the populations of the dots can be 
analytically written as 



61 (t) =e- 2r »"*(l + e 2r ^*)/2 



b 2 (t) 



— P — ik z — 2T 3p t 



*>*(-! + e 2r «*'*)/2 
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From Eq. (9), one realizes that there is always 50% 
chance for the two dots to evolve into the state: ||J.) + 
e -ik Q z a means that, for example, the singlet 

[triplet] entangled state can be created if k z = (2m + l)n 
[2to7t] with m being an integer. The advantage is that the 
entangled states can be generated with a remote sense, 
such that one can manipulate/control one qubit without 
affecting another. One should be reminded again such 
an entanglement is generated from the collective decay 
(super-radiance) in one dimension 22 , not from the trivial 
single-mode Rabi coupling, since the decay rate T sp is 
still present in Eq. (9). 

In summary, we have shown that radiative decay of col- 
loidal QD excitons into surface plasmons can be greatly 
enhanced at certain values of the exciton bandgap. The 
enhancement is due to zero-slope in dispersion relation, 
and one has to treat the decay dynamics with a non- 
Markovian way. In addition, an idea of creating remote 
entangled states between two QDs is also proposed and 
can be tested with current technology 23 . 
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